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ABSTRACT 

Using non-linear time series analysis, along with surrogate data analysis, it is shown that the various types 
of long term variability exhibited by the black hole system GRS 1915+105, can be explained in terms of a 
deterministic non-linear system with some inherent stochastic noise. Evidence is provided for a non-linear 
limit cycle origin of one of the low frequency QPO detected in the source, while some other types of variability 
could be due to an underlying low dimensional chaotic system. These results imply that the partial differential 
equations which govern the magneto-hydrodynamic flow of the inner accretion disk, can be approximated by 
a small number ( « 3 - 5) of non-linear but ordinary differential equations. While this analysis does not reveal 
the exact nature of these approximate equations, they may be obtained in the future, after results of magneto- 
hydrodynamic simulation of realistic accretion disks become available. 

Subject headings: accretion, accretion disks - black hole physics - X-rays: binaries - X-rays: individual (GRS 
1915+105) 
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1. INTRODUCTION 

Black hole X-ray binaries exhibit variability on a wide 
range of timescales ranging from months to milli-seconds. 
This variability is expected to provide important clues about 
the geometry and structure of these high energy sources. 
Moreover, a detailed study of the phenomena may eventually 
be used to test the relativistic nature of these sources and to 
understand the physics of the accretion process. 

The standard first step to quantify variability is to com- 
pute the power spectrum which is the amplitude squared of 
the Fourier transform of the light curve. The power spec- 
tra is expected to give information about characteristic fre- 
quencies of the system which might show up either as spec- 
tral breaks or as near Gaussian pea ks, i.e Quasi-Periodic Os- 
cillations (OPO) (e.g.lBel loni et a/.l i200U ^msick & Kaaret 
1200 It iRodriguez et al\ l2002ii) . The shape of the power 
spectra, combined with the observed frequency dependent 
time lags between different energy bands, have put con- 
straints on the ra diative mecha nism and geometry of th e emit- 
ting region (e.g. iNowak ef g/.l 1999; Misra 2000; Cuil fT99l 
iPoutanen & Fabian 1 999^. 'Nob ih et a/.i2001.) . 

The accretion disk which powers these sources maybe un- 
dergoing stochastic variations at different length and time 
scales. This is expected, especially since the viscosity mecha- 
nism which drives the accretion process is now believed to be 
caused by turbulence induced by magnetic instabilities. These 
structural variations may then manifest as variability in the 
observed intensity after being convolved with the main ra- 
diative process that cools the disk. Additionally, these local 
stochastic variations could couple resonantly with the global 
disk structure and produce coherent oscillations. In most anal- 
ysis, it is implicitly assumed that response of the disk to these 
fluctuations, either by global resonance and/or radiative pro- 
cesses, is linear in nature. However, there is some evidence 
that the response of the disk may be non-linear The ratio 
of the twin high frequency QPO observed in black hole and 
neutron star systems, can be explained in a model where the 
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global disk resonates nonlinearly with the stochastic fluctu- 
ations ( Abramowicz ef a/. 2003). More compelling j"^*^^^ 
independent evidence has been given by Uttlev et al. (i 2005h 
(see also Timmeref a/. (2000) and Thielef a/. (2001)), who 
argue that the log-normal distribution of the fluxes and the 
linear relationship between RMS and flux imply that the re- 
sponse of the disk is non-linear. They show that an ex- 
ponential response can explain these observations. Earlier 
Mineshige ef a/. (1994), had suggested that the behavior of 
the observed fluctuations could be because the accretion disk 
is in self-organized critical state. In all these models the tem- 
poral behavior of the system is driven by underlying stochas- 
tic variations. Moreover, they address the short (< 10 sec) 
time-scale variability of the systems, although many systems 
also exhibit quite dramatic long time-scale variability. 

Variability of a system may not necessarily be driven by an 
underlying stochastic variation. The system can show com- 
plicated temporal behavior if the governing differential equa- 
tions are non-linear and have unstable steady state solutions. 
In other words, although these systems do not have an ex- 
plicit time dependent term in the equations describing their 
structure, they exhibit sustained time variability. In fact, the 
standard accretion disk theory predicts that the disk is unsta- 
ble when it is radiation pressure dominated and when the vis- 
cous stress scales with the total pressure. Numerical hydro- 
dynamic simulations reveal that under such circumstances, 
the disk would under go large ampl itude oscillations around 
the unstable solution dChen & Taamlil994) . These variations 
occur on a viscous time-scale that may be as large as hun- 
dreds of seconds. Specific dynamic models for the tempo- 
ral behavior of accretion disk, like the Dripping Handrail 
(Young & Scarglei ri996h . have been proposed where the ap- 
parent random behavior is actually deterministic. 

The Galactic micro quasar GRS 1915+105 is a highly vari- 
able black hole system. It shows a wide range of long term 
variability (Chenefa/. 1997; Paulefa/. 1997; Belloni ef fl|] 
1997a) which required Belloni ef a/. (2000) to classify the 
behavior in no less than twelve temporal classes. For some 
classes, low frequency oscillations are observed and the light 
curves resemble those from numerical simulations of unsta- 
ble radiation pressure dominated disks based on the standard 
accretion disk theory. However, the viscosity prescription in 
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the standard disk theory is ad hoc, and while it may represent 
approximately the time averages viscosity, it is not expected 
to reproduce the complex time-varying turbulent induced vis- 
cosity expected in an accretion disk. Moreover, there could 
be different types of instabilities in the disk which could give 
rise to strong variability. Nevertheless, the qualitative simi- 
larities between th e results of the simulation and observations 
iTaam et a/.lll997h . suggests that the long term variability of 
this source is due to deterministic non-linear evolution of an 
unstable disk, rather than a manifestation of an underlying 
stochastic process. 

It is important to determine in a more quantitative and 
general method, whether the long term variability is indeed 
due to deterministic non-linear behavior rather than having 
a stochastic origin. That would imply that the temporal be- 
havior of the system can be described by a small number 
of non-linear ordinary differential equations. In other words, 
the complex non-linear partial differential equations that are 
known to govern the hydrodynamical flow, can be approxi- 
mated by a set of ordinary differential equations and hence 
can be more easily studied and understood. Such an approx- 
imation can be obtained, for example, by studying the non- 
linear time evolution of a dominant linear mode. This tech- 
nique which is briefly described in the next section, was used 
to derive the equations for the well known non-linear Lorenz 
system. A more physical method is to study the time evo- 
lution of spatial ly averaged quan tities as in the one zone ap- 
proximation by Paczvnskil Jl983h for thermo-nuclear flashes 
on compact objects. While the rewards of obtaining such ap- 
proximated equations are far-reaching, these methods are in 
general not straight forward and require good physical intu- 
ition, especially since it is not necessary that they can be ob- 
tained for all systems. Hence, evidence for the existence of 
such an approximation, would give the necessary impetus and 
direction to the study of these complex systems. 

Non-linear time series analysis have been used earlier to an- 
alyze X -ray light curves of astrophysical sources. Lehto et oil 
il993h used the correlation dimension technique to analyze 
EXOSAT light curves of several AGN, and found that one, 
NGC 4051, showed signs of low dimensional chaos. Such 
analysis have been undertaken on noise filtered Tenma satel- 
lite data of Cyg X-1 (Unnoefa/. 1990) and on EXOSAT 
data of Her X- 1 (.Voges ef a/.. 1987,: Norris & Matilskv ,1989.) 
and ASCA data of AGN ARK 567 (iGliozzi et a/.ll2002h . In 
an earlier work, we used a modification of this technique to 
study RXTE observations of GRS 1915H-105 and found ev- 
idence that for four of the twelve temporal classes, the sys- 
tem was_consisteiUwith harboring a low-dimensional attrac- 
tor (iMisra et a/.l2004 . 

However, a positive result for the correlation dimension 
analysis, although highly suggestive, is not conclusive evi- 
dence that the light curve of these systems is due to a low 
dimensional chaotic system. The results need to be counter 
checked using other non-linear techniques and compared with 
surrogate data analysis. Thus, in this work, we analyze the 
different temporal patterns exhibited by GRS 1915H-105 us- 
ing correlation dimension analysis and the singular value de- 
composition technique. For each analysis, we undertake sur- 
rogate data analysis as a counter check. To facilitate better 
physical understanding, we draw analogies between the re- 
sults obtained for GRS 1915H-105 and those for a well known 
deterministic non-linear system. 

It should be emphasized that proving the presence of non- 
linear dynamics using a finite time series with mathematical 



certainty is in the very least difficult and perhaps even not a 
well defined problem. Our motivation here, is not to make 
model independent and quantitative non-linear diagnosis of 
the light curve, which as mentioned above is difficult and am- 
biguous. The goal here is to find out if there are features in 
the light curve which indicate (but may not necessarily rig- 
orously prove) that the temporal behavior of the black hole 
system GRS 1915H-105, can be described as a set of ordinary 
non-linear differential equations. Obtaining such equations 
from first principles would be a major break through in our 
understanding of accretion disks around black holes and thus 
any indication that such equations may exist would be useful 
in formulating them. 

In the next section, the well known non-linear system, 
"Lorenz", is introduced and some of the different types of 
non-linear behavior that the system can exhibit are described. 
In §3, some standard non-linear time series analysis is de- 
scribed and as an example applied to the Lorenz system. In 
§4, these techniques are applied to the observed light curves 
of GRS 1915H-105 and in §5 the results are summarized and 
discussed. 

2. THE LORENZ SYSTEM 

The Lorenz model ( iLorenz was developed from 

the Navier-Stokes hydrodynamic equations for the Rayleigh- 
Benard flow, which describes the two-dimensional convection 
of an incompressible fluid in a cell which has a higher tem- 
perature at bottom and a lower temperature at top. The es- 
sential idea is to choose a dominant linear mode that satisfies 
the boundary condition, and substitute this mode back into 
the hydrodynamic equations to obtain the temporal evolution 
of the mode as a set of ordinary non-linear differential equa- 
tions. The choice of the dominant mode is rather arbitrary 
and hence makes this procedure possible only by physical in- 
tuition and/or prior knowledge of the solution from experi- 
ments or numerical simulations. In spite of the physical (and 
in particular) hydrodynamical origin of the Lorenz model, it is 
prudent to be careful about drawing analogies between it and 
other complex hydrodynamical flows like accretion. For this 
work, it is sufficient to state that the Lorenz model is a math- 
ematical set of non-linear equations which have been derived 
from non-linear partial differential equations, using approxi- 
mations which may or may not be valid. These equations turn 
out to be 

Y=RX-Y-XZ 

Z={S/3)Z-XY (1) 

There are three fixed points, or steady state val- 
ues, of the system: {Xu,Yo,Zo) = (0,0,0) and 
{±y/[8/3][R-l],±y/[8/3][R-l],[R - 1]). is a con- 

trol parameter whose value governs the system's behavior 
which is obtained by solving numerically these differential 
equations. For 1 < R < 14, the fixed point at the origin is 
unstable, while the other two are stable. The system evolves 
towards one of the stable fixed points depending on the initial 
conditions. The temporal behavior is transitory in nature. For 
R > 14, all three fixed points are unstable. For large values 
of R > 100, the trajectory along X and F is a closed loop, 
called the limit cycle (Left panel of Figure[3. For any initial 
condition the system evolves towards the limit cycle. For 
intermediate values of R the behavior of the system is more 
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Fig. 1 . — The Dynamics of the Lorenz system. Left Panel is for R = 500 
showing a limit cycle while the right panel is for R = 28, when the system 
exhibits chaos. 
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Fig. 2. — Top panel: The time series X{t) for the Lorenz system with 
R = 28. Bottom Panel: The corresponding power spectrum. 



complex. The right panel of Figure Q shows the trajectory of 
X and Y for R = 28, where the behavior is termed chaotic, 
i.e. infinitesimally close trajectories diverge exponentially in 
time. Thus, the Lorenz system can exhibit different types of 
non-linear behavior 

Often, in a natural system, only one of the several param- 
eters can be detected and hence the complete dynamic infor- 
mation is not directly available. For example, if only the X 
variable of the Lorenz system was known, then from the time 
series for /? = 28 (top panel of Figure|2ji and the correspond- 
ing power spectrum (bottom panel of Figure|5J, the rich com- 
plexity of the Lorenz system is not easily identifiable. How- 
ever, as discussed in the next section, there are non linear tech- 
niques by which the complete (qualitative if not quantitative) 
dynamic nature of such systems can be reconstructed. 

3. NON-LINEAR TIME SERIES ANALYSIS 




X (I) X (I) 

Fig. 3. — The reconstructed dynamics of the Lorenz system using the 
delay embedding scheme. Left Panel is for R = 500 while the right panel is 
for R = 28. 

One of the standard methods to reconstruct the dynamics of 
a non-linear system from a time series, is the Delay Embed- 
ding Technique (Grassberger & Procaccia 1983). Since the 
dimension of the system ( i.e. the number of governing equa- 
tions or variables) is not a priori known, one has to construct 
the dynamics for different dimensions. In this technique, vec- 
tors of length are created from the time series, X{t), by 
using a delay time r, i.e. 

yi= {X{1),X{2)...X{N)) 

V2= iX{l+T),X{2 + T)...X{N+T)) 

VM= X{1 + {M-1)t),...X{N+{M-1)t) (2) 

where M is the chosen dimension. Typically, the delay time 
T is suitably chosen such that the vectors are not correlated 
i.e. when the auto-correlation function of X{t) approaches 
zero. For the Lorenz system, such a technique can effectively 
reconstruct the dynamics as shown in Figure|3] where the X{t) 
is plotted versus X(r + T). Note the similarity between Figures 
|3land[l] with the identification of X(f + r) with F. 

One of the difficulties of this technique is the choice of 
the appropriate delay time, r, which in principle depends on 
the specific nature of the system. To avoid this ambiguity, a 
modified technique has been proposed by Broomhead & Kin^ 
( 1986) where, a matrix constructed of vectors with unit delay 
time, is decomposed into eigenvectors which then represent 
the dynamics. Application of this Singular Value Decomposi- 
tion (SVD) technique to the Lorenz system is shown in the top 
two panels of Figure0]which is again qualitatively similar to 
Figure[n To generate Figure|4] the time series, X{t) has been 
converted to a uniform deviate, which is convenient when the 
analysis is undertaken on data obtained from natural systems 
(see next section). 

To confirm that the dynamical picture obtained by the 
Delay Embedding technique is indeed due to an underlying 
deterministic non-hnear system and not of stochastic origin, 
it is customary to c ompare the results obtained with those 
of surrogate data jTheiler et al I IT992I) . Here surrogate 
data are generated which have practically the same power 
spectrum and distribution as the original time series, but are 
of stochastic origin. The essential idea is to formulate a null 
hypothesis that the data has been created by a stationary 
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Fig. 4. — Singular Value Decomposition (SVD) technique used to recon- 
struct the Lorenz system dynamics immX{t). (a) Using real data fori? = 500, 
(b) using real data for R = 28, (c) using surrogate data for R = 500 and (d) 
Using surrogate data for R = 28. 



Fig. 5. — Singular Value Decomposition (SVD) technique used to recon- 
struct the Lorenz system dynamics from X(t) with noise contamination at 
20% level, (a) Using data for R = 500 with white noise, (b) using data for 
R = 500 with red noise, (c) using data for R = 28 with white noise and (d) 
using data for R = 28 with red noise. 



linear stochastic process, and then to attempt to reject this 
hypothesis by comparing resuhs for the data with appropriate 
reaHzations of surrogate data. iSchreiber & Schnuta (119961) 
have proposed a scheme, known as Iterative Amplitude 
Adjusted Fourier Transform (lAAFT), which generates 
surrogate data that are more consistent in representing the 
null hypothesis (Kuaiumtzis 1999). In this work, the soft- 
ware which uses tha above technique to generate surrogates 
(jhttp://www.mpipks-dresden.mpg.de/tiseanAriSEAN_2.1/index 
has been used. The bottom panel of Figure^] shows the SVD 
decomposition for surrogate data of the Lorenz system. As 
expected, the dynamical phase picture is qualitatively differ- 
ent for the surrogate as compared to the original time series. 
It should be emphasized that a qualitative or quantitative 
difference between a particular kind of surrogate and original 
data, although highly suggestive, does not necessarily imply 
that the original time series has a non stochastic origin. 
Ideally, the null hypothesis should be tested using several 
different kinds of surrogate data, before a concrete conclusion 
can be reached. 

Data from a natural deterministic non-linear system may 
have contamination from stochastic noise, either inherently 
or due to the detection process. Such contamination gener- 
ally makes it difficult to identify the non-linear behavior of 
the system. Apart from the level of contamination, the effect 
on the analysis also depends on the nature of the noise. At the 
same level of contamination (i.e. at the same percentage of 
noise addition) the effect of white noise ( i.e. when the power 
spectrum of the noise signal is a constant) is more than that of 



red noise ( i.e. when the power spectrum, P{f ) °= f~^)- This 
is illustrated in Figure|5J where SVD plots of the Lorenz sys- 
tem with 20% red and white noise addition are shown. While 
the effect of noise on the limit cycle is to broaden the original 
one dimensional curve, its effect on the complex chaotic be- 
havior can be more pronounced. Note that a 20% white noise 
contamination does not allow the dynamics to be effectively 
reconstructed for the Lorenz system with R = 2E. 
html)'^'- often desirable to have a quantitative measure apart 
iiom a qualitative picture of the reconstructed dynamics. Such 
a quantitative measure is the computation of the correlation 
dimension. Briefly, the computation procedure is to choose a 
large number of points in the reconstructed dynamics as cen- 
ters. The correlation function is the number of points which 
are within a distance R from a center, averaged over all the 
centers, and may be formally written as 



Cm(«) 



N Nc 



; lim lim 



\Xi- 



(3) 



where xj are the reconstructed vectors (Eqn 2), is the num- 
ber of vectors, A^^ is the number of centers and H is the Heav- 
iside step function. The fractional or correlation dimension 
D2{M) is defined as 



D.(M).lim^^i^^ 
^ ^ fi^o d\og{R) 



(4) 



and is essentially the scaling index of the variation of Cm{R) 
with R. With these definitions, for a finite duration light curve. 
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Fig. 6. — D2{M) values (filled circles) for the Lorenz system exhibiting 
chaotic behavior ( R = 28). The dashed lines represent D^iM) curves for 
thirty realizations of surrogate data. The D2 (M) saturates to D™' Rs 2 and is 
different from the curves for surrogate data. 



Fig . 7 . — D2 (M) values for the Lorenz system and corresponding surrogate 
with different levels of white and red noise contamination, (a) 20% white 
noise, (b) 50% white noise, (c) 20% red noise and (d) 50% red noise. 



only an approximate value of D2 (M) can be computed. More- 
over for small values of R, Cm{R) is of order unity and hence 
will not show the intrinsic scaling. Similarly, for large values 
of R, Cm{R) will saturate to the total number of data points. 
Thus for a finite data set, an appropriate scaling region in 
R needs to be chosen which does not suffer from the above 
edge effects, to obtain an approximate value of D2(M). Here 
we use the non-subjective algorithm proposed by Misra et 
(|^04) which establishes a region in R which suffers least 
from edge effect errors and computes D2 (M) with an error es- 
timation. A numerical code which implements this algorithm 
is available at http://www.iucaa.ernet.in/'"rmisra/NLD 

For a chaotic system, DolM) « constant = for M 
greater than a certain dimension Mcr- For the Lorenz sys- 
tem exhibiting chaotic behavior (for e.g. when R — 28), 
l^rott & Rowlands (2001) have computed using a large num- 
ber (A^ > 100000) points that Df = 2.05 ± 0.1 and M„ 3. 
Figure |6l shows the D2{M) values for the Lorenz system (for 
R = 2% using 5000 points). As expected D2{M) saturates to 
D™' « 2 for M > Mcr ~ 3. For pure uncorrected stochastic 
white noise D2 (M) = M for all M and no such saturation ex- 
ists. However, a saturated value of D2{M) does not necessar- 
ily imply that the time series is due to deterministic non-linear 
behavior. In fact, colored noise (for which the power spec- 
trum P{f ) °^ /"") exhibits saturation which c an be c omputed 
to beZJf' = 2/(a- 1) (Osborne & Provenza lell989h . Thus it 
is necessary to compare D2{M) obtained from the time series 
with corresponding surrogate data. This is done in Figure |6l 
where as expected the D2{M) for the Lorenz system is signif- 
icantly different from the values for the surrogate data. For a 
deterministic non-linear system the computed Dj"' is a quanti- 
tative measure of the dynamics and the critical dimension Ma- 
is a measure of the number of equations ( for the Lorenz sys- 
tem Mcr = 3) that are required to describe the behavior Con- 
tamination by noise typically increases the computed D2 (M) 
values and makes it more difficult to distinguish the time se- 
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Fig. 8. — Top: The power spectrum and the flux ( count rate) distribution 
for the X-ray light curve of GRS 1915+105 in the p state. Bottom: The power 
spectrum and the flux ( count rate) distribution for corresponding suiTogate 
data. 



ries behavior from surrogate data. This is demonstrated in 
Figure0 where D2{M) values are plotted for the Lorenz sys- 
tem with different kinds and level of noise contamination. 

4. THE BLACK HOLE SYSTEM GRS 1915+105 

Based on RXTE observations, ' Belloni et ai\ (1200 Ih clas- 
sified the various temporal behavior of GRS 1915+105 into 
twelve classes. Here we take a representative data for each 
class and generate continuous light curves in the energy range 
0-30keV binned at 1 sec intervals. The details of the specific 
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Fig. 9. — Singular Value Decomposition (SVD) technique used to recon- 
struct the dynamics using the X-ray light curve for GRS 1915+105 in the p 
state. Left Panel is for real data while the right panel is for a surrogate data. 



data used in the following analysis are tabulated in'Misr a et al.\ 
(|2004). Fi gureOshows the power spectrum of the X-ray light 
curve of GRS 1915-1-105 when the source is in the p state. The 
presence of a low frequency oscillation with several harmon- 
ics is revealed. The flux (or the count rate) distribution during 
this state has a non Gaussian shape. Although the presence of 
harmonics of the oscillation and the flux distribution, suggests 
that the system is non-linear (e.g. Uttlev et al. 2005), they do 
not necessarily imply that this behavior is due to a determinis- 
tic non-linear system. The bottom panel of Figure|S]shows the 
power spectrum and flux distribution for corresponding surro- 
gate data and these are practically indistinguishable from the 
real case. Thus non-linear time series analysis as described in 
the previous section is required to ascertain whether the tem- 
poral behavior is due to deterministic non-linear behavior. 

The SVD reconstruction of the dynamics for the data is 
shown in Figure |5] The presence of a main cyclic feature 
is seen, which is similar to that of the Lorenz system with 
R = 500 (Figure O. This indicates that the observed oscil- 
lation is of a limit cycle origin. Note that the surrogate data 
which has the same power spectrum and distribution has a 
different SVD reconstruction (Figure|9}. 

The power spectra and flux distribution for two other class 
of GRS 1915H-105, 7 and x> are shown in Figure [Tol A 
low frequency QPO with harmonics similar to that seen in 
the p class is detected for the 7 class. During the x class 
no such QPO is detected. However, for both these classes, 
the SVD reconstruction of the dynamics does not show any 
qualitative features (Figure [TTl. Such space filling SVD re- 
constructions are characteristics of systems with uncorrelated 



stochastic noise. In agreement with this, the D2{M) values are 
« M and show no saturation, for these two states (Figure fm i. 
Moreover, the D2 (M) are consistent with those obtained from 
corresponding surrogate data. Data from the (p and S classes 
show similar behavior and are not shown. These results imply 
that the variability of these classes are of stochastic nature, in 
contrast to that of the p class. In particular, the observed low 
frequency QPO in the p and x class have a different origin. 

The behavior of the rest of the classes is more complex. 
For example, the power spectra for A, 9, k and f3 class do 
not show any strong low frequency periodicity (Figure \12i . 
while the flux distribution is typically non-Gaussian (Figure 
\13i . The SVD reconstruction of these data (Figure O show 
structure which qualitatively maybe similar to that of a chaotic 
system with noise (Figure|5} and are somewhat different than 
the SVD reconstruction of the corresponding surrogate data 
("Figure [Tsl. More evidence is provided by the D2{M) values 
for these data sets (Figure [T6l which show saturation and are 
different from the values computed from corresponding sur- 
rogate data. 

5. SUMMARY AND DISCUSSION 

Using non-linear time series analysis techniques, evidence 
has been provided that the black hole system GRS 1915-1-105 
exhibits stochastic variability when it is in x> 7> and S 
classes. In the p class, the low frequency QPO detected in 
the power spectrum can be attributed to a non-linear limit cy- 
cle behavior. For the other classes, the system shows behav- 
ior similar to a non-linear chaotic system with some inherent 
stochastic noise. 

Before considering the implications, it is prudent to em- 
phasize that the results presented here need to be confirmed 
by more comprehensive non-linear data analysis. Since the 
results are based on a null hypothesis test involving surrogate 
data, a more comprehensive test should take into account pos- 
sible different kinds of surrogate data. 

The safiiration of D2(M) to Df « 3-4 (FigureO in anal- 
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Fig. 11. — Left Panels: Singular Value Decomposition (SVD) tech- 
nique used to reconstract the dynamics using the X-ray light curve for GRS 
1915+105 in the 7 and x classes. Right panels: D2{M) values for the 7 and 
X classes compared with corresponding surrogates. 
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Fig. 12.— The power spectra for the X-ray light curve of GRS 1915+105 

in the A, 6, n and f) classes. 
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Fig. 14. — Singular Value Decomposition (SVD) technique used to recon- 
struct the dynamics using the X-ray light curve for GRS 1915+105 in the A, 
8, K and /3 classes. 



ogy with that of the Lorenz system with noise (Figure^), im- 
pHes that the temporal behavior of the system is governed by 
3-4 underlying global equations. On the other hand, it is be- 
lieved that the temporal evolution of an accretion disk can be 



represented by local magneto-hydrodynamic equations which 
are non-linear differential equations in both time and disk ra- 
dius. Thus, an encouraging aspect of the results obtained 
here is that these complicated magneto-hydrodynamic equa- 
tions may be reducible to (or at least approximated by) a set 
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Fig. 15. — Singular Value Decomposition (SVD) technique used to re- 
construct the dynamics using surrogate data corresponding to the X-ray light 
curves for GRS 1915-1-105 in the A, 0, k and f3 classes. 
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of simple (although still non-linear) global equations, which 
are functions of global parameters like size of the disk, aver- 
age density, pressure etc. Such a set of global equations, will 
provide a complete picture of the global properties of the in- 
ner accretion disk and in principle may enable the testing of 
general relativity in the strong field limit. However, obtain- 
ing these equations from the magneto-hydrodynamic equa- 
tions directly is challenging. Perhaps spatial averaging of 
these equati ons to obtain a one (or two zone) models like as 
done by Paczvnski ( 1983) may be the most promising option. 
Such approximations may be obtained with insights gained 
from the results of future sophisticated numerical magneto- 
hydrodynamic simulations of realistic accretion disks. 

Non-linear time series analysis can also be used to compare 
results of numerical simulations with real data. If the sim- 
ulation correctly encompasses the basic non-linear behavior 
of the system, the resultant simulated light curve should have 
similar D2 (M) values and SVD reconstruction. 

In analogy with the Lorenz system, the different temporal 
behavior (corresponding to the different classes) can be at- 
tributed to variation of a single control parameter, which for 
the Lorenz system is R (Eqn[0. For GRS 1915H-105, this con- 
trolling parameter could be the time averaged accretion rate, 
which is determined by the outer boundary condition. One 
of the key attributes of a deterministic non-linear system is 
that the time variability need not be due to a time varying 
parameter For example, the equations defining the Lorenz 
system (Eqn0 do not have an explicit time dependent term. 
For GRS 1915H-105, it is possible that the different temporal 
classes arise from variations in the accretion rate, while the 
variability for a given class is due to an inherent non-linearity. 
Thus it is not necessary to postulate other varying parameters, 
apart from the accretion rate (like magnetic field etc) to ex- 
plain the source's complex variability. If this is true, then the 
average bolometric luminosity which should be proportional 
to the average accretion rate, should be different for each 
class. However, the bolometric luminosity can be estimated 
only by model dependent spectral fitting with the caveat that 
the time-averaged spectra may not be adequately represented 
by steady ones. Nevertheless, the results presented here pro- 
vide the possibility of identifying a single control parameter 
which determines the temporal class exhibited by the system. 

The temporal behavior of GRS 1915H-105 on time-scales of 
hours and months is unique, but may be related to the state 
changes that occur for other black hole systems (like Cygnus 
X-1) on time-scales of months. This difference between GRS 
1915H-105 and the other black holes has led to the speculation 
that GRS 1915H-105 has a rapidly spinning black hole. Identi- 
fying the governing non-linear equations for GRS 1915H-105 
will not only provide insight into the temporal variability of 
other black hole systems, but will perhaps give an explanation 
for its uniqueness. 



Fig. 16. — D2{M) values for the A, 6, k and /3 classes compared with 
corresponding sun'ogates. 
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